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ABSTRACT 

Context. The galaxy cluster XMMU J2235. 3-2557 (hereafter XMM2235), spectroscopically confirmed at z = 1.39, is one of the 
most distant X-ray selected galaxy clusters. It has been at the center of a multi-wavelength observing campaign with ground and 
space facilities. 

Aims. We characterize the galaxy populations of passive members, the fhermodynamical properties and metal abundance of the hot 
gas, and the total mass of the system using imaging data with HST/ACS (1775 and zsso bands) and VLT/ISAAC (J and K$ bands), 
extensive spectroscopic data obtained with VLT/FORS2, and deep (196 ks) Chandra observations. 

Methods. Chandra data allow temperature and metallicity to be measured with good accuracy and the X-ray surface brightness profile 
to be traced out to 1' (or 500 kpc), thus allowing the mass to be reliably estimated. Out of a total sample of 34 spectroscopically 
confirmed cluster members, we selected 16 passive galaxies (without detectable [Oil]) within the central 2' (or 1 Mpc) with ACS 
coverage, and inferred star formation histories for subsamples of galaxies inside and outside the core by modeling their spectro- 
photometric data with spectral synthesis models. 

Results. Chandra data show a regular elongated morphology, closely resembling the distribution of core galaxies, with a significant 
cool core. We measure a global X-ray temperature of kT = 8.6+j 3 , keV (68% confidence), which we find to be robust against several 
systematics involved in the X-ray spectral analysis. By detecting the rest frame 6.7 keV Iron K line in the Chandra spectrum, we 
measure a metallicity Z = 0.26+{J- 2 ° Z Q . In the likely hypothesis of hydrostatic equilibrium, we obtain a total mass of M tot (< 1 Mpc) = 
(5.9 ± 1.3) x 10 14 M G . By modeling both the composite spectral energy distributions and spectra of the passive galaxies in and outside 
the core, we find a strong mean age radial gradient. Core galaxies, with stellar masses in excess of 10" M , appear to have formed 
at an earlier epoch with a relatively short star formation phase (z = 5 - 6), whereas passive galaxies outside the core show spectral 
signatures suggesting a prolonged star formation phase to redshifts as low as z ~ 2. 

Conclusions. Overall, our analysis implies that XMM2235 is the hottest and most massive bona-fide cluster discovered to date at 
z > 1, with a baryonic content, both its galaxy population and intracluster gas, in a significantly advanced evolutionary stage at 1/3 of 
the current age of the Universe. 

Key words, galaxies: clusters: individual: XMMU J2235. 3-2557 - galaxies: evolution - X-rays: galaxies: clusters 



1. Introduction 

Over the past two decades, considerable effort has been devoted 
to discovering ever more dist ant galaxy clusters using differ- 
ent observational methods (e.g. lRosati et alj|2002l for a review). 

* Based on observations carried out using the Advanced Camera 
for Surveys at the Hubble Space Telescope under Program ID 10698; 
the Very Large Telescope at the ESO Paranal Observatory un- 
der Program IDs 274.A-5024(B), 077.A-0177(A,B), 074.A-0023(A), 
077.A-0110(A,B) 



These studies have been traditionally motivated by cosmolo gical 
applic ations of the cluster abundance at high redshift (e.g. lVoid 
2005, for a review), and also by the use of clusters as labora- 
tories to investigate galaxy evolution. Clusters provide a conve- 
nient and efficient way of studying large populations of early- 
type galaxies, which provide stringent tests on galaxy e volution 
mode ls in the current hierarchical formation paradigm dRenzinil 
2006), because they are the most massive galaxies with the old- 
est stellar populations (at least out to z ~ 2). Clearly, the higher 
the redshift the stronger the leverage on theoretical models. In 
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addition, galaxy properties in clusters can be contrasted with 
those in field surveys, which have multiplied in recent years, thus 
extending the baseline over which environmental effects can be 
studied. 

X-ray selection of clusters has been central in these stud- 
ies, as it naturally provides gravitationally bound systems (as 
opposed to simple overdensities of galaxies) with a relatively 
simple selection function. Using ROSAT serendipitous surveys, 
supported by near IR imaging and spectroscopy with 8-10m 
class telescopes, the redshift envelope was pushed to z — 1.3, 
with only 5 cluster s discovered at z > 1 , approximately one 
per square degree (Ro sati et alj 120021) . The extension of the 
same technique to XMM-Newton serendipitous surveys has led 
to the discovery of two cl usters at z > 1.3 to date, XMMU 
J2235.3-2557 at z = 1.39 (iMullis et al.l l 2005h ("hereafter M05 1 
andXM MXCS J2215 9-17 38 atz = 1 .46 dStanford et al.ll2006h . 
See also Lam er etaT] ((2008) for a newly discovered, high X-ray 
luminosity massive cluster at z ~ 1 . 

With the advent of the Spitzer observatory, an alternative and 
efficient way to unveil distant clusters over large areas (as red 
galaxy overdensities in the IRAC and optical bands) has been 
developed. This technique has given notable results, with three 
clusters spectroscopical ly confirmed at z > 1 .3 in the IRAC 
Shallow Survey (ISCS; lEisenhardt et al.ll2008l) and one in the 
SpARCS survey dWilson et al.ll2009l) . 

Detailed investigations of galaxy populations in the X-ray 
luminous clusters at z > 1 have been conducted with HST/ACS 
in combination with the VLT and the Keck telescopes. The 
study of clusters at z = 1.10, 1.24, 1 26, 1.27 dMei et al.ll2006l 
iDemarco et alj 120071 : iMei et all 120091) and th e aforementioned 
systems at z = 1.39 dLidman etall 120081) and z = 1.46 
dHilton et alj 120091) . have revealed tight red sequences for the 
early-type galaxies, with scatters only marginally larger than 
those in local clusters, implying that most of their stellar mass 
was assembled at z > 3 with passive evolution thereafter. While 
a change in the morphology-density relation of early type galax- 
ies (E+S0 galaxies) has been observed at z ~ 1, elliptical galax- 
ies still dominate the cluster galaxy p opulation up to z ~ 1.2 
dPostman et alJl2005l:lHo"lden et al.ll2007l) . 

A comparison of cluster and field early-types of sim- 
ilar stellar mass has revealed a mild but significant dif- 
ference between th e star formation h istories in the two 
environments (e.g. iGobatetalJ 120081: iRettura et all 120081: 
Ivan Dokkum & van der Marell 120071) . a result which is pre- 
dicted by current hierarchical galaxy formation models (e.g. 
iMenci et ai1l2008l) . To date, such comparative studies can only 
be carried out at z < 1 .4, whereas the existence of a substantial 
population of old, massive, passively evolving early -type galax- 
ies in the field is now we ll established up to z ~ 2 dKriek et alj 
120061 ICimattfeFaIll2008h . 

By pushing cluster studies to higher redshifts, where evo- 
lutionary time scales become comparable to the age of the 
Universe at these redshifts, one would expect to detect signif- 
icant evolutionary effects. However, as discussed in this paper, 
this has not been the case so far, even after probing two-thirds of 
the look-back time, not only for the galaxy populations but also 
for the thermodynamical properties and chemical enrichment of 
the hot gas measured with follow-up Chandra observations in 
z > 1 clusters. 

Our current understanding is that relations, such as the 
red sequence in the color-magnitude diagram, the morphology- 
density relation, the Lx - Tx relation for the intracluster gas, 
emerge at z < 2. For example, the study of a proto-cluster 
at z — 2.16 identified around a powerful radio galaxy pro- 



vided evidence of a forming red sequence dZirm et alj 120081: 
Koda ma et alj 120071) . which likely takes 1-2 Gyrs to form 
dGobat et al. 2008). It is unfortunate that such a transition in 
the assembly process of galaxy clusters seems to occur in a red- 
shift range where spectroscopic observation are particularly dif- 
ficult. In this spirit, we have carried out a multi-wavelength study 
of one the most distant clusters known, XMMU J2235. 3-2557 
(hereafter XMM2235) at z = 1,39, which was the first dis- 
tant cluster confirmed (M05) as part of the on-going XMM- 
Newton Distant C luster Project (XDCP. iBoehringer et"al]|2005b 
Fassbender 2007). In this paper, we use spectro-photometric ob- 
servations of XMM2235 in the optical/near-IR and X-ray bands 
to characterize its galaxy population (particularly the passive 
spectroscopic members) and the thermodynamic status of the hot 
intracluster gas and to measure its total mass. In section |2j we 
present the VLT and HST data used in this paper, including an 
extensive spectroscopic campaign which yielded 34 confirmed 
cluster members. In section [3] we model the underlying stellar 
populations of passive galaxies and constrain their star formation 
histories. In section|4j we present deep Chandra observations of 
XMM2235, the methods of analysis and the resulting measure- 
ments of its temperature, metallicity and mass. In section|5j we 
discuss the results. 

H = 70 km s" 1 Mpc~', Q,„ = 0.3, Q A = 0.7 are adopted 
throughout this paper. In this cosmology, 1' on the sky corre- 
sponds to 0.5 Mpc at z = 1-39 

2. VLT and HST Observations 

2. 1 . Optical and near IR data 

Following the relatively shallow discovery imaging data ob- 
tained with the VLT/FORS2 in the R and z bands (M05), 
XMM2235 was observed in the J and K s bands with ISAAC 
on the VLT in January 2006 and August 2007. More recently, 
deeper and wider-field near IR observations of XMM2235 were 
o btained with th e HAW K-I instrument on the VLT and presented 
in lLidman et all d2008l) . 

The ISAAC observations used here consist of a mosaic of 
3x3 pointings, covering 7.5 x 7.5', with exposure times of 
30 and 20 minutes each in the J and Ks filters, respectively 
(program ID 274.A-5024(B)), with the exception of the central 
pointing which was observed for 60 and 45 minutes, respectively 
(ID 077A-0177(A)). The ISAAC data were reduced with the 
ESO/MVM softwarefl Photometric zero points (ZP) were de- 
rived by observing on e or two standards ea ch night from the pho- 
tometric catalogue of lPersson et alJd!998l) . All images were laid 
on an astrometric grid with a pixel scale of 0.15"that was based 
on VLT/FORS2 images. In order to derive a uniform ZP across 
the mosaic, we used independent observations obtained with the 
Sofl camera on the ESO-NTT, also reduced with ESO/MVM, 
for which accurate ZPs were derived with several photometric 
standards observed before and after the observations. These stars 
yielded rms scatters of < 0.02 mag in both J and Ks . The 5x5' 
Sofl field includes the central ISAAC field and covers partially 
the other eight tiles, thus allowing the photometric uniformity 
across the mosaic to be controlled and independently quantified. 
The FWHM of all stars in each field was measured, ranging from 
0.4 to 0.7". Gaussian smoothing was used to obtain a uniform 
PSF in J and Ks mosaiced images (with a FWHM of = 0.51" 
in K5 and 0.73" in J). A sample of 62 stars across the mosaic 



1 http://archive.eso.org/cms/eso-data/data-packages/eso-mvm- 
software-package 
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Fig. 1. Color image of XMM2235 obtained from the combination of i, z (HST/ACS) and K s (VLT/ISAAC) filters. Overlaid X-ray 
contours are from a Chandra 196 ks observation (smoothed by a gaussian with cr = 2"; the lowest contour corresponds to 3<x 
above the background). The image is 2' across, corresponding to 1.0/i^J Mpc at z — 1.4, and is centered on RA = 22 h 35 m 20 s .81, 
Dec = -25°57'40". North is up, east is to the left. 



were selected, using size information from the ACS images and 
their (J - K$)ab color, which was restricted to (J - Ks)ab < 1, 
from the SofI images. These stars were used to to compare pho- 
tometry across all fields, using SofI ZPs as fiducial ones. As a 
result, the ZPs in three ISAAC tiles were adjusted by a few % 
to match the SofI photometry, so as to produce mosaics with a 
uniform ZP, whose accuracy was estimated to be 0.03 mag in K 
and 0.04 mag in J as based on the rms scatter from stellar pho- 
tometry with respect to SofI photometry. We also checked that 
these ZPs are consistent with 2MASS photometry using a small 
number of unsaturated stars in the 2MASS catalog. 

XMM2235 was observed with the ACS Wide Field Camera 
on the Hubble Space Telescope on June 27 2005 in the F7755W 
and F850LP passbands (program ID 1069 8), as part of th e 
ACS Intermediate Redshift Cluster Survey dFord et all 12004). 
Subsequent ACS visits of XMM2235 took place in 2006 as part 
of the HST Cluster Supernova Survey (Dawson et al. 2009). 
This work, however, is based on the first visit only, amount- 
ing to an integration of 5060 sec in £775 and 6240 sec in z&so- 
Analysi s of the full HS T data set are presented in accompanying 
papers dJee et al.ll2009l Strazzullo et al. in prep.). The ACS data 
were process ed with the late s t versi on of the "Apsis" pipeline 
described by Blak eslee et al.l (120031) . We used the AB photo- 
metric system calib rated with the ACS/WFC zero points from 
ISirianni et ail d2005l) . 



The ISAAC mosaic images were aligned with the ACS im- 
ages and a four band (£775, Z850, J and Ks) photometric catalog 
was constructed using SExtractor. Magnitudes were first mea- 
sured in apertures of 0.75" radius in the ACS images and 1" in 
the ISAAC images, then corrected to 2 and 4", respectively, to 
take into account the PSF difference between ACS and ISAAC. 
Such aperture corrections were derived from the median growth 
curve of unsaturated stars in the field. Magnitudes were also cor- 
rected for galactic extinction according to Schlegel (1998). The 
adopted corrections are 0.043 mag in (775, 0.032 in z&so, 0.019 in 
J and 0.008 in K s . A color-composite image of XMMU J2235.3- 
2257 is shown in FigQ] with the Chandra X-ray contours over- 
laid in green (see SectJUi. The peak of the X-ray emission is 
within 2" of the brightest cluster galaxy (BCG) and is elongated, 
from the southwest to the northeast, following the distribution of 
the red cluster galaxies in the core and the major axis of the 
BCG. 



2.2. VLT Spectroscopy 

Three separate observing runs were originally devoted to the 
spectroscopic follow-up study of XMM2235 using FORS2 on 
the VLT. Two runs were carried out using the MXU mode, pro- 
gramme ID 074.A-0023(A) in October 2004 (see M05) and ID 
077.A-0177(B) in July 2006, using two slit masks each. More 
spectroscopic observations were carried out in July 2006, as part 
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of a program to search for distant Type la supernovae hosted by 
early type galaxies in galaxy clusters at z > 1 (Dawson et al. 
2009). For the latter, the MOS mode of FORS2, which consists 
of 19 movable slits that can be moved into the focal plane, was 
used. This allowed rapid follow-up of transient events, such as 
Type la Supernova. In the assignment of slits, transients were 
given highest priority, likely cluster members were given second 
priority and other galaxies were given third priority. Two masks 
were used to follow three transients. Total integration times var- 
ied from 2800 seconds in the first mask to 5400 seconds in the 
second. The number of slits that were used in MXU mode was 
typically more than double the number that were used in MOS 
mode, with greater flexibility in placing targets on slits. All ob- 
servations were carried out with the 3001 grism which has a dis- 
persion of 2.6A/pixel, corresponding to a resolution of ~ 10A 
for the adopted slit width of 1". In combination with the OG590 
order blocking filter, this configuration provides a wavelength 
range of 5900 - 10000A. In the 077.A-0177(B) run, the OG590 
filter was not used thus extending the coverage down to ~5500A. 
For observations with the MXU masks, integration times ranged 
from 4 to 8 hours, depending on the faintness of the targets in a 
given mask. 

Targets were selected using colors and magnitudes, however 
not homogeneously throughout the whole observing campaign. 
Galaxies lying on the R — z red sequence, with zab ^ 23, were 
selected first (M05). Then fainter objects with colors consistent 
with late and early type galaxies at the cluster redshift were tar- 
geted in subsequent masks. As more photometric information 
(ACS-/, z, ISAAC-/, K s ) became available after the first run, we 
adopted new criteria, using Kab < 23 and i - K s , J - K s colors 
ranging from the bluest confirmed star-forming member to the 
reddest member of the red sequence (with a 0.3 mag conserva- 
tive margin). This corresponded to 1.3 < (/77s - K s )ab < 3.7. 

Spectra w ere reduced with stan dard techniques using IRAF 
tasks (see e.g. lDemarco et alJd2007l) ). The RVSAO/XCSAO rou- 
tine was used to measure redshifts via cross-correlations with a 
range of template spectra. Redshifts for faint late type galaxies 
were often derived from the [Oil] line alone. Based on signal- 
to-noise and cross-correlation coefficients (R), spectra were fi- 
nally assigned quality flags of A,B,C depending on whether the 
redshifts was considered secure (R > 4), acceptable {R « 3), 
and unknown respectively. The first three observing runs yielded 
129 redshifts of quality A or B, from 168 extracted spectra, in a 
field covering ~ 3 x 3' (or 1.5 x 1.5 Mpc). From this sample, 28 
redshifts (24 A's+4 B's) lying within +2000 km/s of the median 
redshift zcl = 1 -390 were selected as cluster members, twelve of 
which (1 1 A's+1 B's) are star forming galaxies with a detectable 
[OH]^ 3727 A line (EW([OII])< -5A ). Magnitudes of cluster 
members range from zab — 21.6 (the BCG) to the spectroscopic 
limit of zab — 24. 

More recently, a fourth spectroscopic programme (ID 08 1 . A- 
0759D, PI: M.Tanaka), mostly aimed at the outskirts of the clus- 
ter, was carried out in July-August 2008 with the same 3001 
grism using two MXU masks. J - K s colors were used to tar- 
get galaxies at the cluster redshift over an area of 3 Mpc ra- 
dius. These observations yielded 6 extra passive members, two 
of which are within the central 2' and are therefore included in 
the sample used in this paper. 

In Figf2] we show the redshift distribution 30 cluster mem- 
bers, 18 passive and 12 star-forming galaxies. The latter all lie 
at projected distances larger than 250 kpc. We estima ted the ve- 
locity dispersion using the ROSTAT algorithm of iBeers et alj 
( 1990), which yields cr,, = 802^ km/s with formal boostrap er- 
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Fig. 2. Measured redshifts of galaxies in a field of 3' around 
XMM2235 (with a bin of 0.01). The inset shows a blow-up of 
the distribution (grey histogram in the main figure) around the 
cluster redshift, with a bin Az = 0.0025. The top axis gives 
the rest frame velocity centered on the median value z = 1.390. 
The red histogram includes cluster members within +2000 km/s. 
The hatched region denotes star-forming members. The dashed 
curves is a gaussian with the estimated velocity dispersion of 
802 km/s. 

rors. A detailed dynamical analysis of the entire redshift sample 
will be presented in a forthcoming paper. From the redshift cat- 
alog, we defined a sample of 16 passive galaxies, which are cov- 
ered by the ACS and ISAAC data, and whose spectrophotomet- 
ric properties are modeled in the following section to age-date 
their underlying stellar populations. A complete catalog with 
redshifts from the whole spectroscopic campaign, as well as pho- 
tometric and morphological measurements using the enhanced 
HST/ACS observations, HST/NICMOS and VLT/HAWK-I data, 
will be presented in a forthcoming paper (Strazzullo et al., Nunez 
et al. in prep.). We have verified that the i,z,J,K$ photometry 
used in this paper is in very good agreement with the one in the 
final augmented imaging data set. 

3. Star-formation histories of passive members 

The spectro-photometric data described in the previous section 
were used to constrain ages and star formation histories (SFHs) 
of the 16 passive galaxies (i.e., those with no detectable [Oil] 
emission) covered by the ACS and ISAAC dat a . To t his end, 
we used the technique described in iGobat et al] d2008l) . which 
combines both the spectral energy distributions (SED) and spec- 
tra of galaxies to characterize the underlying stellar populations 
with spectral synthesis models. This generally allows stronger 
constraints to be obtained, by mitigating inherent degeneracies 
among stellar popul a tion p arameters. The same technique was 
used in ISantos et al.l (120091) to age date early-type galaxies in a 
cluster at z — 1 . ACS imaging shows that all the galaxies have 
early-type morphologies, however we defer the morphological 
analysis to another paper (Strazzullo et al. in p rep.). 

We used a grid of Bruzual&Charlot (I2003L BC03) r-models, 
characterized by a delayed exponential star-formation rate (SFR) 
0r(O = (t/T 2 ) exp(-//r), where t is the time since the onset of 
star formation and t a characteristic time scale (corresponding 
to the time of the SFR peak). Solar metallicity and a Salpeter 
initial mass function were adopted. Since spectra extend to rest- 
frame wavelengths shorter than 3200A, we used BC03 templates 
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Fig. 3. Composite spectra of the four brightest galaxies in the 
core (R CL < lOOkpc) of XMM2235 (top) and of the spec- 
troscopically passive members (12) in the outskirts (150 < 
Rcl < 1000) kpc (bottom). Best fit BC03 models to the spectro- 
photometric data of each subsample are overploted in red. The 
ffi symbol corresponds to the atmospheric A band. This region is 
not used in the analysis. 
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Fig. 4. Top: composite spectrum of the 16 confirmed passive 
galaxies of XMM2235 with the best fit BC03 model overplot- 
ted in red. Bottom: an average spectrum of early-type galaxies 
at a similar redshift from the GDDS and GMASS spectroscopic 
surveys in the field. 

which include the iPicklesi (Il998l) stellar library in the mid-UV. 
SEDs are constrained by the four photometric bands described 
in Section im J, Kj, 2775, zsso- Stellar masses for the 16 passive 
galaxies are obtained, together with ages, T, and characteristic 
times, r, by fitting the SEDs with the constraint that T is bound 
to be less than 4.22 Gyr, which limits the star formation epoch 
to z < 15. Stellar masses range between 5 x 10 10 and 5 x 10 1 'M 
apart from the BCG The four red galaxies in the core (within 
100 kpc radius) all exceed 10 n M o with the BCG being by far 
the most massive of the four, with M^bcc ~ 9 x lO n M . 

In order to have a more robust estimate of the SFHs of mas- 
sive passive galaxies in XMM2235, we stacked the SEDs and 
spectra of all 16 members and then separately the four in the 
core and the 12 lying in the outskirts at cluster-centric distances 
150 < R C l < 1000 kpc (Figs [3] and I). Note that these coad- 
ded spectra have a similar S/N ratio, despite the difference in 
the number of galaxies in each sample. This is due to the higher 
luminosities of the core members. The composite spectra have 
a S/N ~ 7, making the fit with BC03 models meaningful. As 



Fig. 5. Star formation histories (in arbitrary units) derived from 
the fit of the composite spectro-photometric data with BC03 r- 
models for the sample of passive galaxies in the core and out- 
skirts of XMM2235. The lower x-axis shows cosmic time; the 
upper x-axis marks the corresponding redshift. 



described in Gobat et al.| d2008l) . we reconstruct the SFH for a 
given subsample by first fitting their average SED, which yields 
a lcr confidence region in the [T, t) space. The subgrid of mod- 
els within this region is then used to find the best fit solution 
to the average spectrum. Finally, the SFH which best fits the 
spectro-photometric data is obtained by computing the mean, 
at each epoch, of all models i(/(T, r) lying within the 3<t con- 
fidence region of the x 1 fit to the average composite spectrum 
and average SED (22 models for the outskirts sample and 13 
for the core sample). The result is shown in Fig|5] A stark dif- 
ference in the derived SFHs of galaxies in the core and the 
cluster outskirts is apparent. Galaxies in the core have formed 
their stars at early epochs, with star-formation weighted ages 
of Tsfr = 3.4*2? Gyr, consistent with a single burst of star 
formation at zf = 5.3^. Those lying outside the core have 
Tsfr - 2.0^ g Gyr, corresponding to zf — 2.7+9 - !?, and a fi- 
nal formation redshift (when 99% of the stellar is assembled) of 
Zfin = 2.0^Q 2- We note that the areas of the SFH curves in Figj5] 
are arbitrarily normalized to one, since we are only interested in 
the difference of the SF time scale between the two samples. It is 
clear that the inferred SFH is model dependent and is generally 
not a unique solution. The difference in the underlying stellar 
populations of the two samples is significant, however, and is 
evident from the inspection of their composite spectra and best 
fit BC03 models (Figj3]). When compared with galaxies in the 
core, galaxies in the outskirts show significant post-starforming 
features (deep balmer lines) as well as a smaller 4000A break 
in addition to having a bluer average SED. Galaxies in the out- 
skirts show residual star forming activity about ~ 1 Gyr before 
the epoch of observation. Note that the [Oil] line is not detected 
even in the stacked spectrum of all the passive galaxies, consis- 
tent with no significant star formation at z = 1 .4. 

These fin d ings a re in very good agreement with those of 
iLidman et al.l (12008b . who used the mean color and scatter of 
a high-quality J - Kg color-magnitude diagram of XMM2235, 
based on VLT/HAWK-I measurements, to infer a formation red- 
shift of zf =4-5 for galaxies in the cluster core (Rcl < 90 kpc), 
in contrast with those outside the core which instead exhibit a 
larger scatter and bluer colors. The SFHs in Fig|5] which was de- 
rived from independent spectro-photometric data, elucidate the 
difference between these two populations. This result is consis- 
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tent with a scenario where the central galaxies formed rapidly at 
high redshift and ceased forming stars early while galaxies out- 
side the immediate core continued star formation for at least 1 
Gyr longer, or alternatively formed most of their stars > 1 Gyr 
later. 

The average spectrum of the 16 passive galaxies and the 
best fit BC03 model to the spectro-photometric data are shown 
in Fig0] For comparison, we also show the composite spec- 
trum of field early-type galaxies at a similar redshift and cov- 
ering are similar range o f stellar masses, whi ch was obtained 
by combining 1 1 G DSS dAbraham et al.ll2004l) and 3 GMASS 
dCimatti et alj |2008) spectra in the range 1 .3 < z < 1.5. It is in- 
teresti ng to not e that, similar to wh at was found by G ob at et ail 
d2008l) (see also lRettura et al.l2008l) . who compared the SFHs of 
early-type galaxies in a massive cluster at z — 1 .24 with those in 
the field, the 4000A break in cluster early-type galaxies appears 
to be larger than those in field galaxies, suggesting a higher for- 
mation redshift or a shorter duration of the star forming phase 
for early-type galaxies in massive clusters. iMenci et al.l (2008) 
found this observational scenario to be in good agreement with 
expectations from semi-analytic models of galaxy formation in 
which SFHs of cluster galaxies of a given mass are modulated 
by environmental effects. 



4. Chandra Observations 

XMM2235 was observed with the Chandra ACIS-S detector in 
VFAINT mode in five exposures of 44 ks (Obs ID 6975), 24 ks 
(Obs ID 6976), 80 ks (Obs ID 7367), 33 ks (Obs ID 7368) and 
15 ks (Obs ID 7404). The data were reduced using the Chandra 
CIAO software V4.1 (December 08 release) and related cali- 
bration files , startin g from the level 1 event file as described in 
iRosati et~ail ((2004). The total effective exposure time is 196 ks 
after the application of this reduction procedure. 

In Fig. Q] we show a color composite image combin- 
ing the HST/ACS and VLT/ISAAC observations with overlaid 
Chandra X-ray contours in the soft band. The diffuse X-ray 
emission from the cluster can be traced out to r ^ 1' (2<x above 
the background), corresponding to 0.5 Mpc. The Chandra im- 
age immediately shows that point sources do not significantly 
contaminate the clu ster emission, as it is sometimes the case i n 
distant clusters (e.g. lStanford et al.l200TtlBignamini et al.l2008l) . 

4. 1 . X-ray Spectral analysis 

We performed a global spectral analysis in an elliptical region 
with semi-major and semi-minor axes a = 42" and b = 38" cen- 
tered on the peak of the photon distribution (RA= 22 A 35 m 20 5 .82, 
Dec=-25°57'40.3", J2000), after masking out point sources. In 
this aperture, which maximizes the S/N, we detected approxi- 
mately 1900 net counts in the 0.3-10 keV band. The background 
was obtained from a large annulus around the cluster position, 
after subtraction of point sources. We also avoided regions where 
missing columns significantly decrease the background inten- 
sity. The background photon file was scaled to the source file by 
the ratio of the geometrical area. The response matrices and the 
ancillary response matrices were computed for each exposure 
with the tool mkwarf and mkacisrmf applied to the extracted 

regions. 

The spectra were analyzed with XSPEC vl2.3.1 dArnaudl 

119961) and fit t ed wi th a single temperature MEKAL model 
dLiedahl et al] dl995lV ). where the ratio between the elements 
are fixed to the solar value as in Ander s~& Grevessel (1989, 



hereafter angr). These values for the solar metallicity have re- 
cently been superseded by the new values of lGrevesse & Sauvall 
(1 19981) . who used a 0.676 ti mes lower Fe solar abundance 
(see also lAsplund et al.l 120051) . However, we prefer to report 
Iron abundance in units of angr since most of the literature 
still refers to these old values. Since our metallicity depends 
only on the Fe abundance, updated metallicities can be ob- 
tained simply by rescaling by 1/0.676 the values reported here. 
We modeled the Galactic abs orption with the tool tbabs (see 
IWilms. Allen. & McCrayfeOQOl) . 

Spectral fits were performed over the energy range 0.5-7 
keV. We exclude photons with energy below 0.5 keV in or- 
der to avoid systematic biases in the temperature determina- 
tion due to uncertainties in the ACIS calibration at low ener- 
gies. We checked that the measured temperature does not de- 
pend significantly on the lower bound, E\, of the energy range 
(by changing E\ in the range 0.3-1.2 keV, T changes by ~ 3%). 
We also exclude energies above 7 keV where the signal is dom- 
inated by noise. We used three free parameters in our spectral 
fits: temperature, metallicity and normalization. The local ab- 
sorption was frozen to the Galactic neutral hydrogen column 
density Nh — 1.5 x 10 20 cm' 2 , as obtained from radio data 
( Dick ey & Lockmanll 1990b . and the redshift was set to z — 1 .39, 
as measured from the optical spectroscopy. Spectral fits were 
performed using Cash statistics (as implemented in XSPEC) of 
source plus background photons, which is preferable for low 
signal-to-noise spectra. We also performed the same fits with^ 2 
statistics (with a standard binning with a minimum of 20 photons 
per energy channel in the source plus background spectrum) and 
verified that our best-fit model gives a reduced x 1 ~ 1 -03 for 
107 degrees of freedom. All quoted errors below correspond to 
lcr, or the 68% confidence level for one interesting parameter. 

The Chandra folded spectra of XMM2235, extracted from 
the whole elliptical area and the inner core, are shown in Fig. [6] 
The corresponding extraction regions are indicated on the in- 
set image. A significant Fe K line is visible at E ^ 2.8 keV. 
Confidence contours from the spectral fit in the Z - kT plane 
are shown in Fig. [7] The fit to the global spectrum gives a best 
fit temperature of kT = 8.6^ 2 keV, and a best fit metallicity of 
Z = O.26^Q 2 gZ (lcr error bars). By leaving Galactic absorption 

as a free parameter, we obtained Nh = 1 ^Iq'^ X 10 20 cirT 2 , thus 
consistent with the Galactic value. In this case, the best fit tem- 
perature is consistent (within lcr) with the aforementioned best 
fit value. Interestingly, leaving the redshift free, a four parameter 
fit yields z = 1.37^q 4 , which shows how accurately the redshift 
can be determined from the X-ray data alone due to the detection 
of the iron line. 

The best fit model corresponds to an unabsorbed flux within 
the elliptical aperture of (2.7 + 0.1) x 10~ 14 erg cirT 2 s _1 in the 
0.5-2 keV band. This corresponds to a luminosity of (2.1 +0.1) x 
10 44 erg s _I /i^ 2 in the rest frame 0.5-2 keV band, and a bolo- 
metric luminosity of (8.5 ± 0.4) x I0 44 h^ erg s _1 . One can use 
the best fit double /J-model of the surface brightness profile de- 
scribed below to extrapolate these luminosities at larger radii. 
For example, values need to be multiplied by a factor 1.3 to en- 
circle the flux within r = Ihji Mpc. Thus, the bolometric lu- 
minosity is Lx,bohm - 1-1 x 10 45 /i^ ( 2 erg s _1 . We note that the 
X-ray luminosity and temperature of XMM2235 are consistent 
with the Lx.bohm 7 T re lation determined for hot clusters at z < 1 
fe.g. lRosatie7aIll2002l) . 

We performed several checks in order to assess the robust- 
ness of our results against systematics. By changing background 
regions, e.g. choosing different external circular regions as op- 
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Fig. 6. Folded X-ray spectra (data points) and best fit MEKAL 
model (solid line) from Chandra observations (196 ks) of 
XMM2235 at z spec = 1.39. The top spectrum is extracted from 
the total elliptical aperture with semi-axis a = 42", b - 38"; 
the one on the bottom from the cluster core (8" radius). The in- 
set shows a cutout (120 X 120") of the cluster in the 0.5-2 keV 
band (smoothed with a <x = 2" gaussian) with the two extraction 
regions. A significant redshifted Fe 6.7 keV line is visible. 



posed to the annulus above, the best fit temperature varies by 
only a few %. More recently, a possible calibration problem of 
ACIS-S at high energies has been under discussion, which has 
supposedly been alleviated by the latest calibration files used 
here. It has been claimed that it biases the temperature of hot 
low-z clusters to higher values, an effect that is limited to 5% 
when the energy range is restricted to 5 keV (L. David, private 
communication). We note that by using Chandra calibrations of 
2007 shortly after the observations, we obtain a temperature of 
9.2 keV (still consistent with the current value at the \cr level), 
whereas the best fit value remains stable over the last two CIAO 
software releases. We also note that the bias is reduced, because, 
at z = 1 -4, the exponential Bremsstrahlung cut off is observed 
at significantly lower energies. By repeating our fit in the energy 
range 0.5 - 5 keV, we find that the best fit temperature decreases 
by only 0.04 keV. 

As shown below, the Chandra data reveal a significant ex- 
cess in the center of the X-ray surface brightness (SB) profile 
of XMM2235, which is usually interpreted as the presence of a 
cool core. A spectral analysis of the inner 8"(67 kpc) radius, 
containing only 400 net counts, yields a best-fit temperature 
kT = 6.7+Jg keV, and an Iron abundance of Z = 0.59^ Z 
(dashed contours in Fi^7]i. This suggests the presence of a cool 
core with higher Iron abundance, a result that needs significantly 
more statistics to be confirmed. We also repeated the spectral 
analysis removing this cool core region of 8.0"radius (the dotted 
l-cr contour in Fig. |7J from the total elliptical region and find 
kT = 9.3^;° keV. 

4.2. X-ray spatial analysis and Mass determination 

The distribution of the soft X-ray emission is clearly elongated in 
the NE direction, resembling the distribution of c luster galaxies 
in the core (see Fig.[T]and also lLidman et alj2008l) . Inspection of 
the surface brightness (SB) profile readily shows the presence of 
excess emission in the central ~40 kpc, which indicates the pres- 



ence of a cool core, a cl ear signature of the dynamically relaxed 
state of the cluster (see Sant os et"ail d2008l) for a discussion of 
cool cores in this as well as in other high-z clusters). Here, we 
report on the modeling of the SB profile, enabled by the angular 
resolution and depth of the Chandra observations, which is then 
used to estimate the mass of XMM2235. 

The presence of a cool core is generally revealed by a gas 
density excess in the central regions of clusters, compared to 
a single y8-model profile, so as to guarantee hydrostatic equi- 
librium. This is clearly illustrated in FigJS] By using a pro- 
jected single y6-model plus background, Y.{r,R c ,B) = [1 + 
(r/R c ) 2 ]~ 3 P +l l 2 + B, we obtain a poor fit in the inner region with 
R c = 90.1 + 12.8 kpc, B = 0.61 + 0.05 with a reduced X 2 of 
2.05. A much better fit is obtained with a double /?-model fit, 
T,(r,R cU 8i) + I.(r,R c2 ,B 2 ) + B, which yields R ci = 45.6 + 11.5 
kpc,/?! = 1.0 (fixed), R c2 = 159.2 +11.2 kpc, B 2 = 0.8 + 0.05, 
X 2 r = l.03. 

This information, assuming hydrostatic equilibrium and 
isothermality of the gas, leads to a deprojected gas density pro- 
file and to an estimate of the total mass. We defer to a forthcom- 
ing paper (Rosati et al. in prep.) a detailed discussion of the X- 
ray mass profile of XMM2235, with a full deprojection analysis 
down to the core and a measurement of the gas fraction, as well 
as independent measurements of the total mass at different radii 
from a dynamical analysis of all cluster members and the model- 
ing of a strong lensing system. We note here that the total grav- 
itating mass profile at large radii (r » r c ) is very well approx- 
imated by a single B— model for the gas density profile which 
leads to the well-known expression: M tot (< r) = 3 - 3 J^* r ' j^s, 
where x = r/r c , fi is the mean molecular weight in atomic mass 
unit (= 0.59), G is the gravitational constant, and m p is the proton 

mass. Hence, M tot (< r) = 1.13 x 1Q H M Q B T(kev) rCMpc)^. 
Using the best fit temperature and the single /3-model param- 
eters above, one obtains a total mass of M tot (< 1 hzlMpc) = 
(5.9 + 1.3) x 1O 14 M . Cluster masses are conveniently measured 
within R500, the radius of the sphere within which the cluster 
overdensity is 500 times the critical density of the Universe at 
the cluster redshift. Using the sc aling relation amon g R500, tem- 
perature and redshift derived in lEttori et ail (120041) one obtains 
for XMM2235 R 500 ~ 0.75 h~^ Mpc and therefore M 500 = M(< 
R 500 ) = (4.4 + 1.0) x 1O 14 M . To estimate the cluster virial 
mass, we have to extrapolate our mass measurement to larger 
radii assuming for example a iNavarro. Frenk. &" White! (119961) 
profile with concentration c — 5. This yields M v i r « M200 = 
I.4M500 ~ 6 x 10 14 M Q . Note that assuming a lower concentra- 
tion, e.g. c — 3, the virial mass would only be a few % lower. 
These values for the total mass indicate that XMM2 235 is the 
most massive cluster at z > 1 discovered to date (see Ettor i et all 
2004 for a compilation of the masses of distant clusters derived 
from Chandra data). 

As discussed in iJee et al.l d2009l) . the X-ray mass profile, 
which is based on the assumptions of hydrostatic equilibrium 
and an isothermal gas, is found to be in very good agreement 
with the mass density profile derived from the weak lensing anal- 
ysis of XMM2235 from deep HST/ACS data, which reveals a 
significant shear well beyond the X-ray SB limit of 1' (or 500 
kp c). Specif i cally, u sing the projection formula for the /3-profile 
in IJee et al .1 d2005), a projected X-ray mass within 1 Mpc of 
M totiPro j(< 1 h~*Mpc) = (9.3 + 2.1) x 1O 14 M is found, in ex- 
cellent agreement wi th the one deriv ed from the shear lensing 
map. We refer to the I Jee etal 1 (120091) paper for a discussion on 
how rare the observation of such a massive cluster is, based on 
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Fig. 7. Confidence contours and best fit temperature and iron 
abundance (crosses) of XMM2235 obtained from the spectral fit 
of the Chandra data (Fig|6]i. Solid contours (1, 2, 3cr confidence 
levels for two interesting parameters) refer to the total elliptical 
region, dashed lines (1, 2cr) to the inner core, dotted line (lcr) to 
the total region with the core removed. 
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Fig. 8. Surface brightness profile from the Chandra soft X-ray 
emission of XMM2235 and best fit beta models. The dotted 
curves indicate the two components of double-/? model (solid 
line) as described in the text. Note the excess of emission in the 
central ~ 40 kpc region respect to a single-/? model. The inset 
shows an X-ray color image (bluer sources have harder spectra) 
of the 2 Mpc region over which the fit is performed. 

the survey volume from which the cluster was discovered and 
the cluster mass function in ACDM cosmology. 

5. Discussion and Conclusions 

We presented a combination of spectro-photometric data from 
VLT and HST, as well as deep Chandra observations of the X- 
ray selected cluster XMMU J2235. 3-2557 at z = 1 .39 and used 
them to characterize the galaxy populations of passive members, 



the thermodynamical properties of the hot gas and the total mass 
of the system. 

Despite surface brightness dimming, due to the high redshift 
of the cluster, the Chandra data show extended X-ray emission 
out to ~0.5 Mpc. The X-ray emission has a regular morphology 
and is clearly elongated in the same way as the distribution of red 
passive galaxies in the core as well as the major axis of the BCG. 
An excess of emission in the inner 50 kpc is clearly detected and 
naturally interpreted as a cool core. The spectral analysis of the 
Chandra data reveals that XMM2235 has a global temperature 
of kT = 8.6+J! keV (68% confidence), which we find robust 
against several systematics involved in the X-ray spectral analy- 
sis. If we assume hydrostatic equilibrium, a well justified condi- 
tion given the relaxed appearance and the canonical Lx/Tx ratio 
of the cluster, and an isothermal gas distribution, the X-ray sur- 
face brightness profile yields a total mass at large radii (r » 100 
kpc) of M tot (< r) ~ 5.9 x 10 14 r(Mpc) M . 

Overall, our analysis implies that XMM2235 is the hottest 
and most massive bona-fide cluster discovered to date at z > 1 . 
These findings are corrobor ated by the we ak lensing analysis 
from deep HST/ACS data dJee et alj 120091) which provides a 
mass profile in very good agreement with the X-ray measure- 
ment. The presence of a significant cool core is addition al evi- 
dence of the advanced dynamical state of the cluster (e.g. Fabian 
11994 . 

In characterizing the passive galaxy population of 
XMM2235, we have extended the analysis of the ACS 
Intermedia t e Red s hift Cluster Survey at 0.8 < z < 1.3 
dMei et al.l 120091; iDemarco et ail 120071) to higher redshifts. 
Using a large sample of spectroscopically identified cluster 
members, we inferred the SFHs of passive galaxies in the 
core (Rcl < 100 kpc) and in the outskirts by modeling their 
spectro-photometric data with spectral synthesis models. We 
find a clear contrast in the underlying stellar populations of the 
two samples. Core galaxies, all with photometric stellar masses 
in excess of 1O U M , appear to have formed at an earlier epoch 
with a relatively short star formation phase (zf = 6-4), whereas 
passive galaxies outside the core show spectral signatures 
that suggest a star formation phase that is prolonged to later 
epochs, to redshifts as low as z ~ 2. The latter are presumably 
the infalling population of galaxies that are in the process of 
populating the red sequence. These results are also consistent 
with their mean col or and scatter of the J - K$ red sequence 
dLidman et al.ll2008l) . 

The on-going star formation is confined to the outskirts of 
XMM2235, as galaxy members with a detectable [Oil] emis- 
sion line avoid the inner 250 kpc region. Also in this respect, 
XMM2235 has already reached an advanced evolutionary stage, 
as star formation is suppressed in galaxies well before they reach 
the center, similarly to low redshift massive clusters. 

We note that the environmental age gradient we find in 
XMM 2235, as well as in other z > 1 clusters (e.g. iMei et al.l 
120091) . is significantly steeper than the relatively mild differ- 
ence found between field and cluster e arly-type galaxies of sim- 
ilar stellar masses dGobat et al.l 12008). This shows the impor- 
tance of environmental effects in driving galaxy evolution, in 
keeping wit h expectations fro m hierarchical galaxy formation 
models (e.g. lMenci et al.ll2008l) . It would be interesting to study 
whether the radial age gradient of galaxy populations becomes 
larger with increasing redshifts. This would require homoge- 
neous galaxy selection criteria in large samples of distant clus- 
ters with homogeneous photometric measurements, a task which 
might be challenging with current data. 
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At increasing redshifts, the observational effects of different 
evolutionary histories of galaxies from cluster to cluster should 
become apparent. Mild differences have bee n detected so far in 
the most distant X-ray luminous clusters (e.g. lMei et al.ll2009l in 
the mean color of the color sequence) o f otherwise very homo- 
geneous populations. iHilton et all d2009l) carried out a detailed 
photometric and morphological study of XMMXCS J2215.9- 
1738 at z — 1.46, currently the most distant spectroscopically 
confirmed X-ray luminous cluster. Based on its X-ray tempera- 
ture and luminosity, as well its velocity dispersion, this cluster 
is less massive than XMM2235 and appears to be in a younger 
dynam ical state based on its velocity distribution (IHilton et al.l 

2007) . Another striking difference between the two clusters is 
the lack of a dominant BCG in XMM2215. However, as in 
XMM2235, the galaxy population of XMM2215 is dominated 
by early-type galaxies with ages of ~ 3 Gyr, based on the scatter 
of its red sequence, implying a major episode of star formation 
at epochs zf ~ 3-5. Unfortunately, it is difficult to directly 
compare these results with those we obtained for XMM2235 as 
the mean ages of the early-type galaxies were estimated with 
different methods. A similar comparative study with optically 
selected clusters would be very interesting but is still lacking. In 
general, the study of cluster-to-cluster variance in distant clus- 
ters, leading to different evolutionary histories of cluster galax- 
ies, could provide further tests for galaxy evolution models. A 
robust comparison among different clusters however would re- 
quire a homogeneous observing campaign with the same set of 
instruments and filters so as to minimize systematics in photo- 
metric measurements and model dependent K-corrections. 

Another notable result of our study is the measured metal 
abundance Z = 0.26^ jgZ Q of the intracluster medium (ICM) 
from the detection of the Fe-K line in the Chandra spectrum, 
essentially coming from the cluster core (Rcl ^ 200 kpc where 
most of the Chandra counts are contained). This extends to even 
larger look-back times the evidence that the ICM was already 
significantly enriched in dista nt clusters, as found in a systematic 
study of clusters at z < 1.3 (Bales traet al.ll200~7l : Maughan et alj 

2008) , and therefore that metal enrichment of the ICM was 
mostly complete at early epochs. In light of the derived SFH 
of the core red galaxies this is not un expected. Chemical evolu- 
tion models of elliptical galaxies (e.g. Pipino & Matteucci 2004) 
predict that chemical enrichment of massive elliptical galaxies, 
which experience a short burst of star formation at early epochs, 
occurs on a time scale of 1 Gyr. Galactic winds distribute metals 
in the cluster core region on a crossing time scale (< 1 Gyr). As 
a result, one would expect the metal enrichment of the ICM to be 
completed by z > 2.5 (2 Gyr earlier than the cluster look-back 
time of 9 Gyr), consistently with the derived SFH of the core 
galaxies. 

It remains somewhat surprising that such a massive cluster, 
which was discovered from a relatively small survey area (11 
deg 2 ), is found with a baryonic content in such an advanced evo- 
lutionary stage at an epoch corresponding to 1/3 of the current 
age of the Universe, both in terms of its galaxy population and of 
the hot intracluster gas. Further results on the galaxy populations 
of XMM2235 as well as its mass distribution, taking advantage 
of an augmented multi-wavelength data set, will be presented in 
forthcoming papers. 

Given its large mass, XMM2235 will also likely be an ideal 
target for a new array of Sunayev- Zeldovich (SZ) experime nts 
which are already operational (e.g. lStaniszewski et alj|2009l) or 
are becoming available. These observations will add indepen- 
dent, important information on the physics of the ICM and the 



total mass of XMM2235 and will be useful to calibrate the ex- 
pectations of a new era of SZ cluster surveys. 
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